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I. INTRODUCTION AND CONCLUSIONS 
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■ Ever since its introduction by Schwinger the heat kernel of Laplace- type operators has become a useful tool 
to deal with one loop effective actions in quantum field theory. This due to the fact that the heat kernel provides a 
manifestly gauge invariant regularization of ultraviolet divergences. An additional virtue, is that, unlike the effective 
action, the heat kernel is a one valued functional. The heat kernel can be applied to study spectral densities of Klein- 
Gordon operators and in the proof of index theorems 0, Q , to compute the ^-function Q and the anomalies of Dirac 
' operators Q , to deal with chiral gauge theories @ and models of QCD , to the Casimir effect Q , to compute black 
hole entropies [jj, etc. Exact calculations of the heat kernel at coincident points are available in particular manifolds 
C^) • [ToL ITT| or for configurations subjected to suitable algebraic constraints (of the constant curvature type) IT3. [Til. [l4| . 
In the general case an asymptotic expansion in powers of the proper time, the Seeley-DeWitt expansion [lallq. is 
• available. The coefficients of the expansion have been computed to rather high orders in several setups, including 
curved spaces with and without boundary, and in presence of non Abelian gauge fields and non Abelian scalar fields, 
^-p . using different methods 0, 0, 0, |2fl|, E3, 0] • The heat kernel expansion at finite temperature has been 

discussed in pH |26| . A generalized heat kernel expansion around non c-number mass terms has been introduced in 
[271 12H |29| . The extension to non commutative quantum field theory has been presented in [3(j • 

The standard heat kernel expansion can be regarded as a double expansion in the strength and in the number of 
derivatives of the background fields. It is therefore most suitable for external fields which are both weak and adiabatic, 
i.e., of slow space-time variation. A resummation of this expansion is provided by covariant perturbation theory [3l| . 
Perturbation theory assumes weak but not necessarily adiabatic fields. In this paper we study a different resummation, 
• . namely, a covariant derivative expansion |32t 13.1 134| . The terms of this expansion have a given number of covariant 
derivatives but any number of scalar fields (these fields playing the role a non Abelian local mass term). So the fields 
are assumed to be adiabatic (and the gauge fields weak, to preserve gauge invariance) but the scalar fields may be 
strong. Using the technique of labeled operators, we are able to write in finite form the contributions to the traced 
heat kernel classified by the number of derivatives. Such contributions can also be regarded as generating functions for 
all the heat kernel coefficients of the standard Seeley-DeWitt expansion with a fixed number of derivatives. Explicit 
results are presented for zero, two, four and six derivatives, involving 1, 1, 6 and 52 summands, respectively. These 
results hold for boundaryless flat space-time but non Abelian gauge and scalar fields. Extension to curved space-time 
should also be possible using the symbols method. A covariant derivative expansion to four derivatives has been 
obtained in [35| for the diagonal heat kernel and in [3(| for the effective action, for curved space-time in the case of 
minimal (i.e. Abelian) scalar field and no gauge connection. 

In Section ^ we define the covariant derivative expansion for the trace of the heat kernel and work out the lower 
order terms. Although the calculation could be done from scratch using the method of symbols, a shorter path is 
provided by the method of Chan ^sj^ which was devised for the effective action. In that Section results are presented 
to four derivatives. We show that after expansion of our results in powers of the scalar field the standard coefficients 
are recovered. Section ITTT1 is devoted to explaining the technique of labeled operators [33l 02 HE HH • Such procedure 
allows to develop a calculus to deal with functions of non commuting variables. The non commutative version of the 
ordinary derivative is shown to satisfy the Leibniz and chain rules, as well as to yield a non commutative version of the 
Taylor expansion. In Section HVI the diagonal heat kernel coefficients, within the derivative expansion, are obtained to 
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four derivatives. This is done through functional variation of the traced coefficients. Unlike the standard expansion, 
in the resummed expansion one finds that each diagonal term comes from the traced term of the same order. This 
yields a consistency condition that is verified by our formulas. Section gives the six derivative contribution. It is 
based on previous results for the effective action found in 39]. Here one important issue is that of finding either a 
short or a systematic expression for the result. This is due to the existence of identities among the possible gauge 
invariant structures, coming from integration by parts and Jacobi identities. In that Section we give a relatively short 
expression for the six derivative contribution, which contains 52 different gauge invariant structures. In Section IVII 
we study the problem of finding a standard basis of structures for the derivative expansion of generic gauge invariant 
functionals. The corresponding problem in the context of the standard heat kernel expansion has been treated before 
by Miiller 0, EH • Some subtleties appear for the derivative expansion because the analogous of the Seeley-DeWitt 
coefficients are now functions (of labeled operators). So, for instance, even if the elements of a basis are complete 
and linearly independent it does not directly follow that the coefficients must be unique; it is at least necessary to 
impose permutation symmetry restrictions to the functions which play the role of coefficients of the expansion. In 
this last Section we construct standard basis of gauge invariant structures for functionals with zero, two, four and six 
derivatives, with 1, 1, 6 and 37 elements, respectively. It is worth noticing that, similarly to what happens for the 
standard heat kernel expansion 30], all the results presented here apply directly to non commutative quantum field 
theory. This is particularly clear when such theories are formulated within the quantum phase space approach |42j . 
which only requires to replace the integral over coordinates by a trace on X-space. Indeed, at no place do we use 
special commutation properties for our symbols, except d v ] = which holds in the non commutative case too. 



II. COVARIANT DERIVATIVE EXPANSION OF THE HEAT KERNEL 



Our goal is to obtain a derivative expansion for the heat kernel of the Klein-Gordon operator 

K = Dl+X. (2.1) 

Here X(x) is a multiplicative operator (i.e. an ordinary function) which is a matrix in some internal space, = 
dfi + V^x) is the covariant derivative, V^x) being a matrix in internal space. Space-time is Euclidean and flat, 
without boundaries and has dimension d. K acts on matter fields in the fundamental representation. 
The standard heat kernel expansion is of the form 

1 °° 

{Ae rK \x) = -^^Y^a n {x). (2.2) 

This is an asymptotic expansion where the a n , known as (diagonal) Seeley-DeWitt coefficients, are r- independent local 
gauge covariant polynomials of dimension 2n. They are constructed with X and in gauge covariant combinations. 
One can choose to order the expansion by powers of r or, equivalently, by the mass dimension carried by the external 
fields Vfj, and X and their derivatives. The quantities d^, and X have dimension 1,1,2 respectively. In spaces with 
boundary, n may take half-integer values |43j, but in our case the index n is a non negative integer. The lowest order 
terms are 

a Q = 1, 
a x = X, 

«2 = \x 2 + + . (2.3) 

In the derivative expansion the terms are classified by the number of covariant derivatives they carry, rather than 
the mass dimension they carry (as in the standard heat kernel expansion) so 

^ TK \ x ) = j^mzZj nA ^), ( 2 - 4 ) 

where A n {x) depends on r and contains 2n covariant derivatives. In this counting, X counts as zeroth order, [D^,X] 
as first order, Z^ v = [D^, D v ] as second order, and so on. (Counting the dimension carried by the background fields 
and by r, A n has dimension 2n.) Technically, the covariant derivative expansion can be defined by introducing a 
bookkeeping parameter A by means of X(x) — > X(Xx) and V^(x) — > AV^(Ax), and then expanding the functional 
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(A 1 x\e rK ( A ) | A 1 x) in powers of A. The derivative expansion is a resummation of the standard expansion, namely, if 
a* denotes the pieces of a n with exactly 2q covariant derivatives (and so with n — q X's) 

A q {x) = Y,r n - q al{x). (2.5) 

n > q 

In what follows we will set r = 1, i.e. remove r from the formulas, since it can be restated at any moment by 
standard dimensional counting. Thus, for instance 

oo 

<*|e*l*>= (4^375 E ^)- ( 2 - 6 ) 
In each A n , X appears to all orders. The prescription to restore r is simply to make the replacement X — » tX, plus 

A n ► T n A n . 

For many purposes it is often sufficient to work with the functional trace of the heat kernel, 

Tre K = J d d xtr(x\e K \x) . (2.7) 

The symbol Tr refers to the full trace on space-time and internal spaces. We will use tr to denote the trace in the 
internal space only. Introducing the short-hand notation 1 

( >-(^/^tr( ) (2.8) 

the standard and derivative expansions take the form 

oo oo 

Tre* = 5>„(:r)> = £<A n (*)> . (2.9) 

n=0 ri=0 

In Tre^ one can use simpler coefficients, b n and B n , which coincide with a n and A n , respectively, modulo by parts 
integration and the trace cyclic property, so that 

(b n ) = (On) , (Bn) = (A n ) . (2.10) 

In particular, 



Tre^ = J2(B n ) . (2.11) 



n=0 

Whereas the coefficients A n and the functionals (B n ) are uniquely defined, there is an ambiguity in the choice of B n 
which is exploited to choose them as simple as possible. The coefficients A n can be obtained from the B n (see Section 
EI). 

The calculation of the coefficients B n can be done using the method of symbols (actually, this methods 
provides A n from which B n is obtained) . In addition to the derivation of the coefficients an important part of the 
calculation is to find a simple expression for them, that is, removing redundancies coming from the trace cyclic property, 
integration by parts and Bianchi identities. Because much work has been devoted to the covariant derivative expansion 
of the effective action, we have found it practical to start from that functional where much of the simplification work 
has already been done. The most useful results for our purposes are found in the work of Chan [32| who obtains 
the derivative expansion of the effective action with a minimum of terms up to four derivatives. This work was later 
extended to six derivatives in |39|. There we find (making explicit terms up to second order) 

Tr l ° gK = J (" l ° SN+ + ' ' ') ' (2 ' 12) 



1 Our notational conventions are summarized in Appendix 1X1 
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where 

N=(k 2 -X)-\ N^ = [D^N]. (2.13) 

In what follows we will consistently use the convention 2 Y^j = [D^Yj] where Yj is an object (such as N, X or Z) 
with an ordered set of Lorentz indices /. So for instance 3 

X MV = [D M , [D v , X]] , Z^ x = [£>„, Z vX ] = [D„, D x }} . (2.14) 

The formula l|2.12(l holds modulo a counter term action which must be a local polynomial (in X, and d^) of 
degree at most d. Such counter terms depend on the renormalization prescription chosen. The derivative expansion 
of the effective action can be written as 

vvtK- vrn^w. <2, 5 , 

J ( 27r ) d ^ T(n + d/2) 

Each F n is gauge covariant, contains a number 2n of D M 's and a number 2n of N's, and has no explicit dependence 
on the space-time dimension d. To four derivatives 



F = -logTV 
1 

2 Afl 



Fi = t;N 2 u , (2.16) 



F 2 = N*Nl -^(N^Nvf - (NN^) 2 -2NZ I1V NN II N V -^(Z^N^ 2 



These terms were obtained in |32j 4 . The six derivative term F% is given in (|5.1fl . It was obtained in |39j| and contains 
45 terms. In this Section we concentrate on terms up to four derivatives and defer the treatment of the six derivative 
terms to Section We will not consider terms with eight or more derivatives in this work. 
First of all we will translate the expansion (|2.15|l to an expansion for the heat kernel: 

= / S"^= / Se4log(*-,) = - / ^log(X-z). (2.17) 



j T 2ni z — K J T 2ni dz J r 2iri 

Here T is a positively oriented path in the z complex plane enclosing the eigenvalues of K. Because the large 
eigenvalues of K lie on the real negative axis, the path is taken starting and ending at — oo. Next we apply Chan's 
formula (|2.15() to Trlog(i^ — z), i.e., with the replacement X — > X — z. The first thing to note is that the counter 
term ambiguities do not survive in the heat kernel, since a polynomial in z does not give a contribution to the integral 
in (|2.17|) . Taking Tr in l|2.17|l and inserting (|2.15|) , one finds that all non explicit dependence on fc M comes in the form 
X — k 2 — z. Making the shift z — » z — k 2 removes all dependence on fc M in F n . (Such change of variables is justified 
at the level of asymptotic expansions we are considering.) Straightforward momentum integration gives then 



that is, 



B n = - I ^e z F n (2.19) 
■ 2m 



r 



(perhaps modulo integration by parts and cyclic property). In this formula the quantities F n are given by the same 
expressions (|2.16|l where now 

N^(z-X)- 1 . (2.20) 



2 Here and else wh ere in this work Y stands for a generic matrix-valued function. 

3 Note that in loll the convention Yin = [D^, Yj] is used ins tea d. Also F M „ there corresponds to %Z^ lv here. 

4 Note that the sign of the fourth term of F2 is incorrect in |32l . 
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The integration over z is easily done for the zeroth order term (undoing the steps in Q2.17fl ) 

Bo = I ^e*lo g N= [ ^^—= e *. (2.21) 
J r 2m J r 2m z — X 

For Bi, we first expand the covariant derivative using the identity 

= NX^N (2.22) 

so that z appears only in N outside covariant derivatives 

F 1 = ^N 2 X^N 2 X^ (2.23) 

(exploiting the cyclic property to move the last N to the first place). 

When this Fi is inserted in (|2.19|l the integral over z cannot be readily done because z appears in two places and 
the operators do not commute in general. Here we apply the technique of labeling the operators [33l 134 . l37l l38| : 
relative to the two X^ in (|2.23|) there are three positions, namely, the position 1 at the left of the first (leftmost) X^, 
the position 2, in between the two X^ and the position 3, after the second X^. (Operators at position 3 can be moved 
to position 1 by the cyclic property.) The operators relative to which the positions are defined (the two X^ in this 
case) are named "fixed operators". The other operators are then labeled according to their position relative to the 
fixed operators and moved to the left (or to any convenient location in the expression). In this way we can rewrite 
(l2~2l|) as 

Fl = \n 2 N 2 X^ = \ J -J— J -J— Xl . (2.24) 

The labels 1 and 2 indicate at which position the labeled operators Ni and N 2 (or X\, X2) should be inserted (relative 
to the fixed operators X^X^). The point of following this procedure is that the labeled operators are effectively c- 
numbers: similar to the time-ordered product or the normal product, they can be written in any order since their 
true position is given by their label. Because they are c-numbers, nothing prevents us from doing the z integration 
as for ordinary functions. We refer to Section ITTT1 for details on the use and properties of labeled operators. Their use 
is crucial for the rest of the paper. After integration over z we obtain 

(B 1 ) = (f(X 1 ,X 2 )Xl) (2.25) 

where X\ is X located at position 1 and X2 is X at position 2. f(Xi,X 2 ) is the ordinary function f(x,y) evaluated 
at X\ and X 2 , with 



1 f dz 1 1 e x — e v 1 e x + e v 

f{x,y) = --^—e*j—^ ^^^=— . ( 2.2(i ! 



(By definition T always encloses the poles, in this case at x and y.) 

The right hand side of (|2.25|) can be regarded as the generating function for all the Seeley-DeWitt coefficients with 
two covariant derivatives. Indeed, making a series expansion in X one obtains (using the cyclic property in the second 
equality) 



Lh • = ( ( -i - i(* + x 2) - ± { xl + xl) - L Xx x* + --)x* 



Xl XXl - —X 2 Xl - —XX U XX U + ■■■) . (2.27) 

12 12 40 60 I 1 ' 

To the order shown, this reproduces the pieces with two covariant derivatives in 63, 64 and 65 17] . 
For the four derivative term we proceed similarly. Expanding the covariant derivatives of N one finds 



F 2 = (N 2 X^N 2 X^) 2 - -{N 2 X^N 2 X v f - A(N 3 X X ^) 2 - (N 3 X^f 
-4N 3 X^NX^N 3 X UU - 2N 2 X f ,N 2 X u N 2 Z flu - \{N 2 Z^f 



N 2 N 2 N 2 N 2 (X' 2 ) 2 - ^N 2 N 2 N 2 N 2 (XM 2 - 4N 3 N 2 N 3 N4Xl) 2 - N 3 N 3 X 2 ^ 
-AN 3 N 2 N 3 XlX vy - 2NlN 2 NlX, t X v Z^ - \n 2 N 2 ZI„ . (2.28) 
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The integrals over z involved in the computation of the B n are of the form 

dz 



I ri , r2 ,...,r n (Xi,X 2 ,...,X n ):= / — e *N?N?.--N?. (2.29) 

Jr 2vri 

They can be computed from the basic integrals 

I°{X u X 2 ,...,X n ):= [ ^e*N 1 N 2 ---N n = J2e Xi U^A 1 r (2-30) 
Jr 27r * ^1 .jr. Xi - Xj 



using 



" 1 / g 

Ir 1 ,r 2 ,...,r n (Xi,X 2 , ■ ■ ■ , X n ) = ^ - — (^gY~J In(Xl, X 2 , . . . , X n ) . ( 2 -31) 

(Alternatively, J rii7 . 2i ... irn (Xi, X 2 , . . . ,X n ) can be obtained from Zj? _| |_ r taking the first n arguments to be Xi, 

then the next r2 arguments to be X2 and so on.) Note that the functions I ri ,r 2 ,...,r n are everywhere analytic in the 
n-dimensional complex plane, and invariant under a common permutation of labels rj and arguments Xj. The basic 
functions I® are completely symmetric and they can be obtained from the following recurrence relation 

= e Xl , (2.32) 

r f Y Y Y Y \ — ^n(Xl, X n -i,X n ) - In(Xi, . . . , X n -i,X n+ i) 
J n+ll A l' ■ ■ ■ ' A nj ^~n+l) — T; 7; ■ 

We will use the short-hand notation 

(X 1 ,X 2 ,...,X n ). (2.33) 
With this notation the lowest terms in the derivative expansion are 

(Bo) = (h) , 

(Bi) = (-\l 2 , 2 xfj, (2.34) 

(B 2 ) = ( (-/ 2 ,2,2,a + 4/3,1,3,1) [Xlf + i/2 A 2,2(V„) 2 + ^3,3^ 

+ 4i3 i i,3-X'^X„„ + 1I 2 , 2 2,X jl X v Z lxv + —I 2l2 Z^ v 
Expanding in powers of X one obtains for B 2 

(B 2 ) = 

+ Y8Q X3Z ^ U + Y25^ 2 ^^^ + YgTj {XX ^X V Z ^ v + X^XX l ,Z^ v X tl X v XZ^ v ) + — XyXw 



which reproduces 62, 63, &4, b 5 to four covariant derivatives [l7|. 

The expressions (|2.34ll . together with the similar expression l|5.2|l for B3, are the main result of this work. Note 
that, compared to the standard heat kernel expansion, in the covariant derivative expansion the numerical coefficients 
of the standard expansion are replaced by coefficients which are functions of labeled X's. 

It can be observed that the coefficient functions found are directly consistent with the cyclic property. For instance, 
the identity 

{fiX^X^X^X^XlXl) = (f(X a ,X i ,X 1 ,X 2 )XX) (2.36) 
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shows that one can always choose the coefficient function of (X 2 ) 2 to be invariant under the cyclic permutation 
(Xi, X 2 , X 3 , X 4 ) — > (X s , X4, Xi, X 2 ) and this symmetry is explicit in (B 2 ) given in l|2.34[l using the permutation 
symmetry properties of the functions I ri ,r 2 ,...,r n - 

There is another symmetry also realized in the B n which will play an important role in what follows. This is mirror 
symmetry, that is, the symmetry under transposition defined by linearity plus the rules 

(AB) T = B T A T , D^ = D^, X T = X. (2.37) 
They imply [D^, Y] T = — [D^, Y T ] and thus 

= (-ir^x...^ , c..^ = (-ir- 1 ^..^ • (2-38) 

In practice, because the total number of Lorentz indices is always even, an equivalent rule is to pick up a minus sign 
for each Zf il ...a n in the expression. 

Mirror symmetry holds for the Klein-Gordon operator K and for the heat kernel e , and is manifest in Bq, B\ and 
B 2 . E.g. 

(f(X 1 ,X 2 ,X 3 )Xf l X vu ) T = (f(X 4 ,X 3 ,X 2 )X uv Xl) = (f(X 3 ,X 2 ,X 1 )X^X vv ) (2.39) 

shows that the coefficient function of X^X VV (namely, 13.1,3 in this case) can be chosen even under transposition of 
Xi and A3. It is curious that the coefficient function of X^X^Z^y has a greater symmetry than required by cyclic 
and mirror symmetries. 

III. LABELED OPERATORS 

In this section we describe useful properties of labeled operators. In an expression with labeled operators, there are 
fixed operators relative to which the positions are defined and labeled operators which carry position labels, e.g. 5 

f(Ai,B 2 ,C 3 , . . .)XY ■ ■ ■ (3.1) 

X, Y, etc, are the fixed operators in this case, A is to be inserted before (to the left of) X (position 1), B between X 
and Y (position 2), C just after Y (position 3) and so on. Such an expression can be defined in two equivalent ways. 
First by writing / as a sum of separable functions, e.g. 

f{x\i x 2 , x 3) . . .) S c nit n 2 ^ n3i ,,Xi x 2 2^3 , 

f(A 1 ,B 2 ,C 3 ,...)XY--- := Yl <>m, na , na ,...A ni XB n *YC n *-- . (3.2) 

7li,7l2,n3,... 

Another example would be the expansion of / as a combination of plane waves, through its Fourier transform. 

Alternatively the expression in (|3.1I) can be interpreted through its matrix elements. Taking \n,A) as a basis of 
eigenvectors of A with eigenvalue a„, and similarly for B and C, and being (n, A\, etc, the corresponding dual basis, 

(n,A\f(A u B 2 ,C 3 )XY\r,C) = Y,f( a ^ m ,c r ){n,A\X\m,B){m,B\Y\r,C). (3.3) 

m 

The important point is that the expression in l|3.1f) depends only on the function / and not on how it is expanded. 

The usefulness of the labeled operators stems from the fact that they are effectively c-numbers since e.g. A\B 2 = 
B 2 Ai, and so they can be used in several ways. For instance, [A, ) can be written as Ai—A 2 , since [A, X] = [Ay— A 2 )X, 
then f([A, ]) can be represented as f{A\ — A 2 ). The well-known identity 

e^ A - 'A = e A Xe~ A (3.4) 

becomes trivial using labeled operators 

e [A, ] X = e A t -A 2X = e A l£ -A 2X = e A Xe -A (3 5) 



5 The symbols A, B, X, etc, are generic and do not refer to those of the heat kernel throughout this Section. 
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This kind of properties have been used in [.3.1 138| to easily derive commutator expansions. 

Labeled operators appear naturally in non Abelian expansions. If one needs to compute 6 f(A + B) to first order 
in B, where A and B do not commute, a standard technique is to transform this problem into e A+B by means of a 
functional transform, and then apply Dyson's formula 



e A + [ dse sA Be^ A + 0(B 2 ). 
Jo 



e A+B = e A + ds e sA Be^~ s > A + 0{B 2 ) . (3.6) 



Using labeled operators one can go further and obtain 



5 '' ; = « *- / 1 dse sAl+(1 -^ A2 B + 0(B 2 ) = e A + ^ ~^j~ B + 0(B 2 ) (3.7) 







(B being the fixed operator). Undoing now the functional transform yields the useful relation, for a generic function 
/(*), 

f(A + B) = f(A) + f{Al J ~ f } A2) B + 0(B 2 ) . (3.8) 

Ai — A 2 

This relation does not rely on the exponential function and can also be established by using, e.g., a Taylor series 
expansion of f(x). 

From the previous expansion we learn that under a first order variation 5 A of the operator A (the function / being 
unchanged) 

5(f(A)) = Vf(Ax,A 2 )SA (3.9) 

where the operation V is defined by 

V/( gl ,x 2 ) := -ft**)- /(*») . (3 . 10) 
xx - x 2 

V maps a one-argument function f(x) into a two-argument function Vf(x,y). Note that S can be any variation of 
A, including e.g. a derivative 

d fl f(A) = Wf(Ax,A 2 )d tl A (3.11) 

or a commutator, S = [X, ], 

[XJ(A)] = Vf(Ax,A 2 )[X,A] (3.12) 

(mdeedVf(Ax,A 2 )[X,A}=Vf(Ax,A 2 )(A 2 -Ax)X = (f(A 2 )-f(Ax))X = [X,f(A)]). 

The operator V generalizes the ordinary derivative to the non Abelian case. Of course, when A and 5 A commute the 
right-hand side of l|3.9|) becomes f'(A)SA (applying L'Hopital rule). It is straightforward to verify that the operator 

V satisfies a Leibniz rule 

V(fg)( Xl ,x 2 ) = Wf{ Xl ,x 2 )g(x 2 ) + f{ Xl )Wg(xx,x 2 ) , (3.13) 

where f(x), g(x) are possibly non commuting (matrix valued) functions of a single variable. Furthermore, V also 
complies with the chain rule 

V(f o g)(xx,x 2 ) = Vf(g(xx),g(x 2 ))Vg(xx,x 2 ) (3.14) 

(in this case f{x) may be matrix valued) as is readily verified. 
In the case of several variables one may need partial derivatives, e.g. 

6(f(Ax,B 2 )X) - Vxf(Ax, A 2> B 3 )6AX + V 2 f(A 1 ,B 2 , B 3 )X5B + f(Ax,B 2 )6X , (3.15) 



By /(A), and similar expressions, we mean an ordinary function f(x) evaluated at x = A, A being an operator, in the analytical 
extension sense, as in e A (and not to completely general operator- valued functions of A). More generally, / may be matrix- valued 
taking values in a different space than that of A. 
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where V i indicates that it acts on the £-th argument of /. 

It is also convenient to define a V operator acting on the space of functions of any number of variables, in such a 
way that if maps a n-variable function f(x\, . . . , x n ) to a (n + Invariable function, as 



V:=^V fe (3.16) 



k 

(Vfc acting on the fc-th argument of the function) and so 

t— t r i \ f{ x li ■ ■ ■ > x k+l, ■ ■ ■ , %n+l) — f(Xl, ■ ■ ■ , %k: ■ ■ ■ , x n+l) /017 > 

V/(xi,...,ar n+ i) = > , 3.17) 

f— Xk~ Xk+l 

where xt indicates that the £-th argument is missing. The definition is such that 

6 (f(A u . . . , = Vf(A u A 2 , . . . , A n+1 )(6A) n (3.18) 
for 6(A) = 8 A and 8(8 A) — 0. Using this operator we can write down the non Abelian version of Taylor's formula 

f(A + B)=e BV *f(A) (3.19) 
(Va emphasizes that V acts on the ^-dependence of the expression). That is 

f(A + B) = f(A) + Vf(A 1 ,A 2 )B + ^V 2 f(A u A 2 ,A 3 )B 2 + ■ ■ ■ 

Ai — A 2 



f(Ai) , f(A 2 ) f(A 3 ) 



(A 2 -A 1 )(A 3 -A 1 ) (A 1 -A 2 )(A 3 -A 2 ) (A x - A 3 )(A 2 - A 3 ) 



B 2 + ---. (3.20) 



In the Abelian case the coefficients reduce to those of the standard Taylor expansion. 

To finish this Section we note a very important point when using labeled operators, namely, that of the regularity 
of the functions of the type f(A\, A 2 , ...) at the coincidence limit of two or more arguments. These functions must be 
regular (free from poles) at the coincidence limit to define meaningful operators. For instance, an expression of the 
type 

' B (3.21) 



A x -A 



is only formal as it refers to any solution Y of the equation [A, Y] = B. Depending on the operators A and B such 
solution either does not exist or is not unique. The operation V (cf. eq. (|3.1U|l ) does not introduce singularities 
(poles) in that limit and the functions appearing in the expansion l|3.20|l are all regular. 

IV. DIAGONAL HEAT KERNEL COEFFICIENTS 

The diagonal coefficients A n (x) can be computed using e.g. the method of symbols, however, having the B n it is 
simpler to derive them from the relation |(| 

S x Tre D l +x = Tr (e D l+ x 8X) (4.1) 

where Sx is a first order variation with respect to X. That is 

ww-^S- (4 - 2) 



For the Seeley-DeWitt coefficients this implies 
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whereas for the derivative expansion coefficients gives 



(4.4) 



This relation allows to obtain A n from B n , but B n can also be obtained from A n though H2.10(l . This implies the 
consistency condition 



It can be shown that this consistency condition is satisfied by our expressions 

/ri,r 2 ,...,r„ = e Xl f ' (X2 — X\, . . . , X n — X\) 



(4.5) 

This only requires the property 

(4.6) 

(this relation codifies the heat kernel property e K — ► e a e K under the shift X — > X + a, a being a c-number). The 
consistency condition is not sufficient to determine B n since it will be satisfied too by all heat kernel-like operators of 
the form exp(X + /(£>£)). 

Using the results in Section ITTT1 relative to the manipulation of labeled operators, one can carry out the functional 
derivative with respect to X indicated in l|4.4|) for (B n ) given in (|2.34|) . However, in the present case it is simpler to 
go back to F n (i.e., before integration over z and labeling of operators) and do the variation there, using the identity 



S X N = NSXN . 



(4.7) 



We illustrate the method with A\ 

1 



Sx(Fi 



(NSX^X^Xf, + N 2 SXNXf,N 2 X^ + N 2 (6X)^N 2 X fl ) 

= ((N 2 X^N 2 X^N + NX^N 2 X^N 2 - (N 2 X^N 2 )^ 5X) 
= ((-2N 2 X„NX fi N 2 - N 2 X^N 2 )8X) . 



(4.8) 



I.e. 



S{Fi 



1 



(2N 2 N 2 N 2 X^X ll + N 2 N 2 X^ 



SX (47r) d / 2 

Following the same procedure with F2 and carrying out the integration over z, we find 



(4.9) 



A 
Ax 
A 2 



2/2,1,2 ZftvZfiv 



+ (4/ 2 , 



2,2 



8/; 



3,0,3 



+ 8/ : 



2, 1,3 J ^fifiu 



x v 



16/3 



+(4/2,1,2,2 - 16/3,0,1,3)^^ 
+(32/3,0,1,3 — 16/2,1,1,3 - 8/2,1,2,2 

+2/3,3 X^ 

+(2/2,3,2 + 4/3,1,3 + 4/2,2,3 

+ 8/3,1,3 XpvXpv 

+(16/3,1,3 + 8/3,2,2)^,^ 

+ (4/2,2.2,2 + 16/3,1,1,3 + 16/2,2,1,3) XpX^Xv 
+(2/2,2,2,2 + 8/3,1,1,3 + 8/2,2,1,3) XpXwXfj, 

+ (4/2,2,2,2 + 16/3,1,1,3 + 8/2,2,1,3 + 8/2,3,1,2 + 8/3,1,2,2 + 8/3,2,1,2) X^X ll X v 
+ (32/3,1,1,3 + 16/3,1,2,2)^,^-^1/ 

+ (4/2,2,1,2,2 + 16/3,1,1,1,3 + 16/2,1,2,1,3 + 8/2,2,2,1,2 + 8/2,1,3,1,2 + 16/2,2,1,1,3) X^X^XjjXu 

+ (4/2,2,1,2,2 + 16/3,1,1.1,3 + 16/2,2,1,1,3) X^XyX^Xjj 

+ (4/2,2,1,2,2 + 16/3,1,1,1,3 + 16/2,2,1,1,3) X„X V X V X^ . (4.10) 
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Ai is a direct transcription of (|4.8|l . However, in A2 we have shortened the expression by using a standard convention, 
namely, we have identified every term with its mirror conjugate and have used only one of the two forms. In other 
words, each term stands for the semi sum of itself plus its mirror conjugate. For instance, 

4/2,2,3 X^X VV := (2/2,2,3 + 2/3,2,2) Xfi^Xw , 
16/3,0,2,2 Z^X^Xy := 8/3,0,2,2 Z^X^Xn + 8/2,2,0,3 X^X^Z^v . (4-11) 

In obtaining A2 we have used the Jacobi identity (here Y represents an arbitrary quantity) 

Y flu = Y v ^ + [Z flu ,Y] (4.12) 

to reduce the number of terms, e.g. by eliminating terms of the type X^ V Z^ V (as h[Zu V , X\Z^ V ) and by canonically 
ordering the Lorentz indices. 



V. SIX DERIVATIVE TERMS, B a 

The term F3 in the expansion (|2.15l) has been computed in |39|: 



F3 = + ^-N a N ct NpNpN J N 1 - 2N a NpN a N 1 NpN 1 + ^ N a Np A 7 N a Np A 7 + 2N a N a NpN 1 N 1 Np 

-ANaNaNpNjNpNy - WNN a NpN a N^N^ + 16NN a NpNpN^N ai + 6N 2 N aa NpNpN 71 
-16NN aa NNppN^N 7 - 8NN a N a pNN 7 N y p - 16NN a N a pNpNN^ - 8NN a NpNN aj Np^ 
+8NN a NpNNp 7 N ai + 10NN a NppNN a N^ - ANN a Npp NN 71 N a + 8NN a NppN a NN 71 

-^NN aa NNppNN 17 + 12N 2 N aa NpN Npjy + ZN 2 N a ppN 2 N ail 

-8Z a pNjNN a NpNN-y + ANZ a pNN 1 N a NpN 1 - 8NZ a pNN a NpN 1 N 1 - 8NZ a pNN a N 1 NpN 1 
—ANZ a pNN a NjNryNp - 16NZ a pNN a NN^Np - 16Z a pN 2 N Ty NN a Np - \QZ a pN 2 N 1 N ia NNp 
-16NZ a pN 2 N ai N 7 Np - 16N Z a pN 2 NaNjyNp - 8N 2 Z a pN 2 N ai Np 1 - lW 2 Z a pN 2 N a Np 17 
-\QNZ a pN 2 N ail NNp + 8N 2 Z aa pN 2 NpN 11 

-2N 3 Z aa pN 3 Z 7 p 7 + 16Z a pN 3 Z jaj N 2 Np + N Z a pN N^N Z a pN N~ t - ±NZ a pNN a NZp 7 NN 7 
-4N 2 Z a pN 2 N 7 Z a pN y + 8Z a pN 2 N 7 Z aj N 2 Np + \&Z a pN 3 Z^N^NNp - 2NZ a pN 2 Z a pNN y N 7 
-4NZ a pN 2 Z aj NNpN y + ANZ a pN 2 Z aj NN^Np - 8Z a pN 3 Z a pN 2 N 11 

~N 2 Z a pN 2 Z a7 N 2 Zp y . (5.1) 

In this formula the number of explicit terms has been reduced by identifying terms related by i) cyclic permutations 
and ii) mirror symmetry. As noted before some conventions here differ from those in 39] . In order to obtain the heat 
kernel coefficient B3 from F3 we use the procedure of Section [n] as in e.g. H2.28fl . This yields 



(B 3 ) — ( — Sl^XappXajy — 24I i i4 : X a Xp a Xpj 1 + (— 12/324 — 12/4i4)A Q A ( g / 3A Q77 
g 

+ -^l333X aa XppX 77 + (— 24/ 411 4 — 24/ 412 3)A Q X a XpXp 11 — 12I i n i X a XpX a Xp 11 

+ (-12/1324 - 6/1414 - 6/2224 + 16/2233 + 4/2323 + 16/ 3 13 3 )^a^a^/3/3^ 77 

+ ( — 24/i4i4 — 24l23\\)X a XpX a pX 11 + 8l3232X a XpX a ~ f Xp J — 8l3232X a XpXp~ ( X a ~ ( 

— 24:Ii4i4X a XpX 7 pX 7a + (— 24/i4i4 — 24/2314 + 16/3223)^0^0/3-^^77 + 8/3232-^0-^0/3-^7-^7/3 

—24:l4i4iX a Xp a X 7 Xpj + (— 6/1414 — 12/1423 — IO/2323 — 8/2332 )X a XppX a A 77 
+ (— 24/13214 — 24/14114 — 24/ 2 22i4 + 32/22313 — 24/ 2 3ii4 + I6/23213 + 32/31223 + 32/31313 
+16/32123)-^o-^o-^/3-^/3-^77 

+ ( — 48/41141 — 48/4i23l)^Q^a^/3^7^/37 + ^l32V32X a X a XpX 1 X 1 p 

+ (— 16/13232 — 16/22232 + 32/31322 — 48/41141 — 48/41231)-^a-^a-^/3-^/37-^7 

+ (— 24/11324 — 24/11414 — I6/21233 — 40/21323 — 24/ 2 i4i3 — 16/31322)^0-^0-^/3-^77-^/3 
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+ (— 247i4ii4 — 24/ 2 3ii4 + 16l32i23)X a XpX a XpX Tr 

+ (16/32222 + 16/32312 — 48/4ii4i)X Q X/3X Q X 7 X ) 3 7 

+ ( — 16/21323 + 10>l32132)X a XpX a X 1 X 1 p + \§l-s23\2X a XpX 1 X a Xp 1 

+ (— 24Jii4ii4 — 48/n4i23 — 8/121323 — 16/l22223 + 32/i23i32 + I6/213123 + 64/213132 

20 64 

— 24/ 2 i4i22 — ^-^222222 + -^"^313131 )X a X a XpXfiX 1 X 1 

+ ( — 48/ii4H4 + 8/121323 + 16/l22223 + 16/l23123 — 16/l32312 — 48/214113 — 16/222213 

+4/222222 + 3'2I-},Y3212)X a X a XfiX 1 XpX 1 
+ ( — 48/n3214 — 24/ii4H4 + 8/i23123 — 32/212313 — 40/213213 — 2l222222)X a X a XpX 1 X 1 Xp 
+(16/i 21323 — 8/213231 + 16/221322 + 2/222222 + 8/312312 — IZIilUl^XaXfjXaXyXisXj 

2 

+ (8/l21323 — q ^222222 )X a XpX 1 X a XfiX 1 



+ 8/323^0^/3/3^707 + (16/233 + 16/323)^0^/377^0/3 + 8/323^0/3^7/3-^07 
— 16/3123^Q^Q^/3-^7/37 + 16/3123^0^/3^70-^/37 

+ 16/ 3 23iX a X a / 3 X 7 Z / 3 7 + (I6/2223 + 32/2313 + I6/3123 + '&21- i 2\3)X a XpX 1 $Z al 

+ (16/l233 + 16/2313 + ^h2l3)X a XpX 11 Z a i3 + (32/3123 + 32/3132 )X a Xfj a X 1 Zf3 1 

+ (16/2223 + 16/2232 + 16/2313 + 16/3123)X a X/3/3X 7 Z Q7 

+ (8/22222 + 32/31123 + 32/31132 + 32/31321 + l&l32l3l)X a X a X{3X 1 Zp- l 

+ (8/22222 + 32/31123 + 32/31132 + I6/31213 + I6/31222 + ^l3213l)X a XpX a X 1 Zpr 1 

+ (32/22123 + 32/22132 + 16/22213 + 4/22222 + 32/23H3 + 32/31123 + 8/3i2l3)X Q ,X ( gX ( gX 7 Z Q , 7 

+ (8/l3231 — 4/22222 + 8/3i2l3)X Q ,X/3X 7 X Q Z/3 7 



+ 2/33^00/3-^7/37 + 16/3l3X Q Z Q /3Z 7 /3 7 + 8/i33X aa Z/3 7 Z/3 7 + (I6/1133 + 2l2222)X a X a Z/3 1 Zp 1 
+^l2222X a XpZ al Zf)~ t + (— 16/1313 — 4/2222)X a X/3-Z/3 7 Z Q7 + (4/2222 + 8l3\3l)X a Z a pX 1 Z fj-f 
+ (4/l331 — -^2222)X Q Z/3 7 X Q ,Z / 3 7 



+ 2^222-^0/3-^07^/37^) ■ (5.2) 

For the proper interpretation of this formula, it is important to recall that mirror symmetric terms have been 
identified, as in (|4.11|) . 

-F3 contains 45 terms but only 40 different structures of fixed operators, since some of the terms differ only by the 
position of underivated iV's. Similarly, B% contains 52 different structures of fixed operators (e.g. X a ZpyX a Zpy), 
each with a coefficient function written as a combination of functions I ritT2 ... r „. Counting each of these as different, 
-B3 contains a total of 147 terms. 

As explained at length in [3Ej , the expression of F3 is not unique due to integration by parts and the Jacobi identity 
(|4.12|) and this is also true for B3. The previous expression for B3 comes directly from F$ using only cyclic and mirror 
symmetries to reduce the number of terms. This is also the shortest expression we have found for B3 from the point 
of view of the number of terms. No systematic minimization of the number of terms in B3 has been attempted (it was 
done in (3^ for F3), nevertheless the existence of a much shorter expression seems unlikely. Alternatively, one can 
try to reduce the number of structures. The length defined from the point of view of the number of structures can be 
reduced from 52 to 37. This is because, as discussed below, all functionals of the type of B3 can be written using a 
standard basis of structures with 37 elements. (Conceivably, a concrete functional such as B3 could be written using 
a smaller number of structures, but this is unlikely.) The price to pay for a smaller number of structures is to increase 
the number of terms from 147 to about two thousand terms. 



VI. BASIS OF STRUCTURES WITH TWO, FOUR AND SIX DERIVATIVES 



Instead of finding shortest expressions, it is also of interest to find a standard basis [401 141| of structures. We discuss 
that problem in this Section. 
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Specifically, one would like to express a generic gauge invariant functional F(D,X) constructed with and X 
and with a fixed number of covariant derivatives, as a combination of structures Tj (the basis, independent of the 
functional F) with F-dependent coefficients FW which are functions of labeled JC's: 

N 

(F)=J2(F {z) (Xi,X 2 ,...)T l ) (6.1) 

i=i 

(the number of arguments in being the number of fixed operators in the structure Ti). 

A standard basis is subject to some requirements, namely, i) the structures in the basis must be sufficient to express 
any functional and ii) all of them must be necessary (i.e., no one can be removed from the basis without spoiling the 
requirement (i)). 

The expressions for Bq, B\ and B 2 in l|2.34[l suggest that the following are standard basis for gauge invariant 
functionals with cyclic and mirror symmetry: For zero derivatives, 

Ti = 1 . (6.2) 

For functionals with two derivatives 

Ti = Xl , (6.3) 

and for four covariant derivatives, 

T\ = {X 2 ) 2 , T 2 = (X^Xv) 2 , T 3 = X 2 ^ , T 4 = X 2 X UU , T 5 = X^X^Z^ , T 6 — Z 2 ^ . (6.4) 

This is indeed so. For instance, for two covariant derivatives, in addition to T\ = X 2 , one could write down a further 
structure X^, however, this is redundant since 

(f(X)X^) = <-[£>„, = {-Vf{X x ,X 2 )Xl) = (f^(X 1 ,X 2 )T 1 ) . (6.5) 

As we have noted at the end of Section lllll the coefficient functions F^ (Xi, X 2 , ...) must be regular in the coincidence 
limit of two or more arguments to define meaningful expressions. If this requirement is not attended one finds that 
formally a smaller number of structures would suffice. For four derivatives these are T\ and T 2 . For instance, T 3 can 
be reduced to T\ as follows 

(f(X 1 ,X 2 )X 2 ^) = (f(X 1 ,X 3 )(X 1 -2X 2 +X 3 )(X 3 -2X 4 + X 1 )D tl D^D v D u ) 
t!Y ^ (X 1 -2X 2 +X 3 )(X 3 -2X 4 + X 1 ) 

f(x lt x a ) {X2 Xi)[X3 X2){Xi _ Xs){Xi —X,X,X V X V 

>( 1 )(X 1 ,A 2 ,A 3 ,A 4 )T 1 ) . (6.6) 

Such reduction is faulty as the identity used in the second step — ► (X 2 — X\)~ 1 X^ is only formal (as noted above, 
the equation [Y, X] — X^ does not imply Y — D^). 

In addition to being sufficient and necessary, for a standard basis one may ask whether the coefficients F^ are 
unique. In general they will not be unique. For instance, for the structure X 2 , to any given F^ (X\, X 2 ), one can 
add an arbitrary odd function f(Xi,X 2 ) = —f{X 2l Xi). Such an addition (f(Xx,X 2 )XT) vanishes identically using 
cyclic symmetry. Therefore, F^ (Xi, X 2 ) can only be unique if one imposes the further condition that it must be 
symmetric under transposition of X\ and X 2 . It is clear that under such restriction F^'(Xi, X 2 ) is unique. (This can 
be verified using the technique of bare structures introduced below.) Equivalently, a functional (F^ (Xl, X 2 )XV) is 
identically zero if and only if the symmetric function F^'(Xi,X 2 ) is identically zero. 

In the four derivative case, the coefficients can always be chosen to have the following symmetries 

77,(1) _ _ _ 

^1234 — ^1432 — ^3214 _ r 3412 ) 

77,(2) _ 77,(2) _ p (2) _ p (2) p (2) _ p (2) _ p (2) _ p (2) 

^1234 — ^1432 — ^2143 ~~ ^2341 — ^3214 — ^3412 — ^4123 ~~ ^4321 ) 

77,(3) _ p (3) 
r 12 — r 21 > 

77,(4) _ p (4) 
^123 — ^321 

17,(5) _ p (5) 
r 123 ~ r 321 > 

F$=f£\ (6.7) 
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where we use a shorthand notation e.g. ^432 = FW(X 1 ,X 4 ,X 3 ,X 2 ). For instance 

F^X^XvZ^u) = — (^F^ZfivXyX^j = — (f^XvX^Z^ = {F^X^XuZ^J^ (6.8) 

(using, in the first equality, that a term and its mirror symmetric have been identified). With these symmetry 
restrictions, these F^ can be shown to be unique. 

We finally come to the six derivative case (always assuming gauge invariant functionals with cyclic and mirror 
symmetries). The 52 structures appearing for B3 in (|5.2|l are neither necessary nor sufficient. A standard basis is as 
follows 

T\ = X a/ 3 7 X a/ 3 7 , T2 = X a Xp^X a fjy , T3 = X a pX a7 X/3y , T4 = XaXaXffyXpry , 

T5 = X a XpX 7 X a {}~ ( , Tq — XaXpXa^Xp-! , Tj = XaXpXp^Xa^ , Tg = X a X a pX 1 X 1 j3 , 

Tg = X a Xp^X a Xf3^ , Tin = X a X a XpX 1 Xp 1 , T\\ = X a XpXpX 1 X ai , T12 = X a XpX a X 1 Xp 1 , 
T13 = X a XjjX 1 X a Xp 1 , Ti4 = X a X a XpXpX 1 X 1 , T15 — X a X a XpX 1 XpX 1 , 
Ti6 = X a X a XpX 1 X 1 Xf3 , T17 = X ct XpX a X 1 XpX 1 , Tig = X a XpX^X a X^X~ ( , 

Tig — XaXpjZpa^ , T20 = X^X^Z^ , T21 = A Q A^X-,,.^^ , T22 = AqX^Xq^Z^ i 

?23 = XaXpXp^Zaj , T24 = XaXafiXjZp^ , T25 = X Q X a X^X 7 Z^ 7 , T26 = X a XpX a X 1 Zp 1 , 

T27 = X ol XpXpX 1 Z oll , T28 = XaXpX^XaZpj , 

T29 = Z a p 1 Z a p 1 , T30 = X Q Z / 3 7 Z Q( g 7 , T31 = X a pZ al Zp 1 , T32 = XaXaZ^Z^ , 

?33 — X a XpZ a jZpy , T34 = XaXpZpyZay , T35 = X a Z Q ^X 7 Z / 3 7 , T36 = X a Zp 1 X a Zp 1 , 

T37 = ZafiZa^Zp^ . (6.9) 

To establish that this set is sufficient we follow the ideas put forward by Miiller in 41] for the standard heat 
kernel expansion. Consider the set of all possible structures with six derivatives, removing those that are redundant 
using i) dummy indices, ii) cyclic symmetry and iii) mirror symmetry. There are 211 such distinct structures. Using 
integration by parts one can always remove all structures where Lorentz indices are contracted within the same factor, 
e.g. the index a in X^X^X^. Next, one can use the Jacobi identity (|4.12|) to choose the order of the covariant 
derivatives within each factor, for instance, if X a ^ 7 X Q( g 7 is retained, X Q , (( 3 7 X Q , 7( g becomes redundant. For the same 
reason all structures of the type (• ■ • K.. Q ^... 7 • • • Z... a p) are also redundant. Finally, XaX^X^Zp^ can be reduced to 
T21 using the Bianchi identity 

^Q/37 — Zp a y + Zy/3 a . (6.10) 

This produces the set of structures in l|6.9|) . 

The above constructive procedure shows that the 37 structures are sufficient. Before showing that they are also 
necessary, we need to discuss the symmetries of their coefficient functions F^\ These are as follows 

_ p(3) _ P (3) _ p (3) p(4) _ p (4) p (5) _ p (5) 

12 — r 21 I r 123 — r 213 ~~ r 321> r 1234 — r 3214' r 1234 — r 4321' 

. p(9) _ p (9) 
' ^2143 ~~ ^3412' 



pW - pW p( 3 J _ p( 3 ) _ p( 3 ) _ pW p( 5 ) _ 

r l2 — r 21 > ^123 — ^213 — r 321' r 1234 — r 3214' r 1234 " 

P (6) _ p (6) p(7) _ p (7) p (8) _ p (8) p (9) _ p . 

^1234 ~~ x 3214> x 1234 ~~ r 3214> x 1234 — ^3412' ^1234 — r i„ , 

_ p(13) _ p (13) p (14) _ p (14) _ p (14) 

^12345 — ^32154' r 12345 ~ r 54321' ^123456 — ^165432 — r 321654' 

P (15) _ p (15) p(16) _ p (16) _ p (16) p (17) _ p (17 

^123456 — r 321654i r 123456 — ^321654 — r 456123' r 123456 — r l 

p (18) _ p (18) _ p (18) p (20) _ p (20) p (27) _ p (27) 

^123456 — r 654321 — r 234561' ^123 — ^321 ) r 12345 — r 54321' 

p (28) _ p (28) p (29) _ p (29) „(31) _ p (31) 

^12345 — r 5432lJ r 12 — r 21 i r 123 ~ 1 213 > r 



_ p(17) _ p (17) 
— ^165432 — ^456123' 



F (28) _ p (28) p (29) _ p (29) „(31) _ p (31) „(32) _ „(32) 

12345 — J 54321' r 12 — r 21 ' ^123 -^213 ' r 1234 — r 3214' 

p (33) _ p (33) p (34) _ p (34) „(35) _ p (35) p (36) _ p (36) _ p (36) 

^1234 — r 3214> r 1234 — r 3214' r 1234 — ^3412' ^1234 — ^4321 — r 3412' 

F (37) _ ^(37) _ ^(37) (6 n) 

r 123 ~ r 213 ~ r 231 ■ I - 11 ; 

We have indicated only the generators of the symmetry group (e.g. for F^ it follows that this function is completely 
symmetric under permutations). As discussed before, only after imposing the symmetries can one expect the coefficient 
functions to be unique for a given functional. A subtlety that did not appear in the two or four derivative cases is 
that one has to take into account not only true symmetries but also quasi- symmetries. For instance, 

7 123^31/ = (f[23 XafiZa-fZfij) = ( F^ 2 ^ Zp 1 Z a ~ f X a p \ = (F 2 il X a pZp^Z a ^ 
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^213* X a pZa^Zfr + (X\ — X2)F^3 Z a fjZ al Zfj 1 



F£l } T 3 i + (X, - X 2 )f£1 } T 37 ) , (6.12) 
(31) 

therefore the antisymmetric component of F^ 23 under transposition of 12 can always be traded by a contribution to 

(31) 

T37 and one can require -F 12 3 to be symmetric. 

To verify that the 37 structures are necessary we have used the following device. We consider a generic expres- 
sion F of the type (|6.1|l with unspecified coefficient functions . Then F is expanded in terms of bare structures, 
namely, structures formed with operators D^, as in the first equality of l|6.6|l . There are five such bare struc- 
tures, DaDaDpDpD-yD^, D a D a DpD 1 D ! 3D 1 , D a D a DfjD 1 D 1 D / 3 1 D a DpD a D 1 DpD 1 , and D a DpD 1 D a DpD 1 (mod- 
ulo cyclic and mirror symmetries) . The reason to do this is that an expression written in terms of bare structures is 
zero if and only if the corresponding coefficient functions vanish (after imposing the appropriate symmetry restrictions 
to those coefficients). That is, there are no identities (like Jacobi or integration by parts) in terms of bare structures, 
so two expressions are equal only if their (symmetrized) coefficient functions are equal. To see that a given structure 
Ti is necessary, i.e., that it cannot always be written in terms of the other structures, it is enough to expand F as a 
power series of X in terms of bare structures and equate it to zero. If Tj were redundant, for any choice of F^ there 
would be choices of the other coefficient functions so that the equation F = would hold true at each order in the 
series expansion. It can be verified that this is not the case for any Ti when one considers the equations at order X . 

In the six derivative case, we have not found a closed proof that the symmetrized coefficient functions associated 
to an expression are really unique. In principle there could exist non trivial identities, that is, sets of non vanishing 
functions F™ producing a vanishing expression F. (If the F™ were not symmetrized or the 7j were not all necessary, 
this would certainly be the case.) To investigate this issue, we have considered large classes of functions F™ of 
the type encountered in B3, i.e., obtained by linear combination of functions I ri ,r 2 ,...,r„, with adjustable numerical 
coefficients. The corresponding expression has been expanded in terms of bare structures and equated to zero. No 
nontrivial identity has been found. Our conjecture is that the symmetrized coefficient functions corresponding to a 
given expression are unambiguous. 
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APPENDIX A: SUMMARY OF THE NOTATION AND CONVENTIONS 

In this Appendix we collect several notational conventions used in the text. 

K = Dl+X, D^ = d IM + V fl , Z^:=[D^D V \. (Al) 

1 

(47rr) d / 

a n -> r n a n , B n (X) -> T n B n (rX) . (A3) 
Ym = [D„ Yj]. (A4) 



< > : =7tAt777 I d d xtr( ). (A2) 



Units restoration: 



Indices convention: 



Symmetric functions: 



(rt - 1)! \dX t/ 
t=i y ' v ' i=l jV 



Mirror transformation: 

AB > B T A T , Xfj, 1 ___f ln > X^...^ , ^ 1 ...^ 1 „ ► ~Zf ll „,f ln ■ 
Mirror symmetry convention: 

Y:=\{Y + Y T ). 
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